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ABSTRACT 

We study the PAH emission from protoplanetary disks. First, we discuss the dependence of the PAH band ratios on 
the hardness of the absorbed photons and the temperature of the stars. We show that the photon energy together 
with a varying degree of the PAH hydrogenation accounts for most of the observed PAH band ratios without the need 
to change the ionization degree of the molecules. We present an accurate treatment of stochastic heated grains in a 
vectorized three dimensional Monte Carlo dust radiative transfer code. The program is verified against results using 
ray tracing techniques. Disk models are presented for T Tauri and Herbig Ae stars. Particular attention is given to 
the photo-dissociation of the molecules. We consider beside PAH destruction also the survival of the molecules by 
vertical mixing within the disk. By applying typical X-ray luminosities the model accounts for the low PAH detection 
probability observed in T Tauri and the high PAH detection statistics found in Herbig Ae disks. Spherical halos above 
the disks are considered. We show that halos reduce the observed PAH band-to-continuum ratios when observed at 
high inclination. Finally, mid-IR images of disks around Herbig Ae disks are presented. We show that they are easier 
to resolve when PAH emission dominate. 
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1. Introduction 

The thermal structure of the surface of protoplanetary 
disks, and in particular the part which gives rise to the mid- 
IR emission, may be strongly influenced by the presence of 
polycyclic aromatic hydrocarbons (PAH). Such molecules 
have small heat capacities and when they absorb individual 
photons their temperature fluctuates strongly. Prominent 
signatures of PAHs are given by detection of their infrared 
emission bands (see Tielens (2008) for a review). Spatially 
resolved mid-IR observations give evidence of gaps and spi- 
ral structures (Jayawardhana et al. 1998; Augereau et al. 
1999). The origin of the PAH emission in protoplanetary 
disks is confirmed for most of the sources (Habart et al. 
2004; Lagage et al. 2006; Doucet et al. 2007; Geers et al. 
2007), although PAH emission is observed at several thou- 
sands of AU away from the star (Siebenmorgen et al., 1998). 
The detection rate of PAH bands in protoplanetary disks 
for T Tauri stars is as low as 10%, for Herbig Ae stars as 
high as 60% using ISOSWS (Acke & van den Ancker, 2004) 
and goes up to 70% using the superb sensitivity of Spitzer 
(Acke et al., 2010). In the disks of the more luminous Herbig 
Be stars the PAH detection rate is below 35% (Acke et al., 
2011). 

The destruction and survival of PAH in the disks of T 
Tauri stars is discussed by Siebenmorgen & Kriigel (2010, 
SK10). The PAH detection statistics in protoplanetary 
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disks seems to be related to the excitation mechanism of 
the PAHs and their destruction, in particular, by hard pho- 
tons (soft X-rays). The PAH destruction is counter-acted 
by the survival of the molecules which might be caused by 
turbulent mixing within the disk. We wish to corroborate 
these claims by detailed radiative transfer models of proto- 
planetary disks around T Tauri and Herbig Ae stars. 

This paper is structured as follows: We first study the 
dependence of the strength of PAH bands on the hydro- 
genation coverage of the molecules and their excitation by 
photons of different hardness hv (Sect. [2]). It is found that 
the spectral shape of the exciting radiation, in particular 
when hard photons are involved, should not be neglected 
when incorporating PAHs in radiative transfer models. A 
method to include stochastic heated grains such as PAH 
in a vectorized three dimensional (3D) Monte Carlo (MC) 
dust radiative transfer model is presented in Sect. [3] As a 
benchmark test we compare results of the MC procedure 
to a ray tracing code. Starlight heated, so called passive 
disks, are treated, which are in hydrostatic equilibrium with 
gas and dust coupled. Disks are irradiated by the stellar 
photosphere and by harder photons with energies of sev- 
eral eV up to 2keV. Vertical mixing of the grains within 
the disks is considered. Flaring disks are computed with a 
one dimensional (ID) radiative transfer code provided by 
Kriigel (2008). Disks are also computed with a MC code 
in 3D, where the combined hydrostatic equilibrium and ra- 
diative transfer equations are solved following the proce- 
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dure by Siebenmorgen & Heymann (2012, SH12). The lat- 
ter method was developed for disks including equilibrium 
heated grains. In those disks ring-like structures and gaps 
appear in the surface structure. In Sect. 0] we present re- 
sults for disks including stochastic heated dust particles. 
The spectral energy distributions (SED) are computed for 
pure disks and disks with dusty halos. The surface bright- 
ness distribution in the PAH bands is discussed. Our main 
findings are summarised in Sect. [5J 



2. PAH excitation and band ratios 

Variations in the relative strength of the 11.3 (im C-H to 
the 6.2 and 7.7/j,m C=C vibrational modes of PAHs are 
collected for a large sample of galactic and extragalactic 
objects by Galliano et al. (2008). The authors conclude 
that the PAH ionisation degree accommodates the observed 
variations of the band strengths (Draine 2011). We point 
out a degeneracy between the ionisation degree and the 
hydrogenation coverage of the molecules in explaining the 
band ratios. We consider a single PAH of constant: ioni- 
sation degree, abundance (25 ppm [C in PAH]/[H]), and 
size. The PAH is made up of 150 C and 37 H atoms. The 
absorption cross section per C atoms in the optical/UV 
is K = 7 x 10~ 18 cm 2 with a cut-off at A cu t = 0.6/im 
(Schutte et al., 1993). We consider 17 emission bands and 
take for simplicity Lorentzian profiles. Parameters are listed 
in Siebenmorgen et al. (2001). They are calibrated using 
ISO spectra of starburst nuclei. The PAH model represent 
an average between neutral and ionised species. Variations 
of the band carriers which may be caused by chemical and 
radiation processing of the molecules in the particular en- 
vironment are not treated. Nevertheless, in disks around 
Herbig stars the centers of the PAH features shift to redder 
wavelengths as the stellar effective temperature decreases 
(Keller et al. 2008, Peeters et al. 2002). In our models such 
a shift can only be compensated ad-hoc. Objects with a 
strong 3.3/^m band are easier fit including an additional 
smaller PAH component or a PAH size distribution. The in- 
fluence of a bimodal size distribution and the impact of the 
hydrogenation degree of PAHs on their mid-IR emission 
is discussed by Siebenmorgen & Kriigel (1992). However, 
models including a single PAH with size as above appears 
sufficient to fit SEDs of various galaxies (Siebenmorgen & 
Kriigel 2007). In order to keep computing times of the MC 
procedure reasonable we treat PAHs of one size. 

First we consider PAH excitation by mono-chromatic 
radiation and study the influence of the hardness of the 
source spectrum on the band ratios. The emitted radia- 
tion is at constant strength (flux F — 100ergs _1 cm~ 2 ), its 
hardness is varied using photon energy between hv = 2 eV 
and 50 eV. As F cx n 7 hv is constant the number of pho- 
tons n 7 decreases as the photon energy hv increases. For 
the same strength of the radiation field the PAH tempera- 
ture distribution function P(T) becomes narrow when ex- 
cited by the many soft photons and widens when excited 
by the few hard photons impinging on the molecule (Fig. 3 
in SK10). 




F(6.2)/F(11.3) 

Fig. 1. Strength of the CC stretching modes at 6.2 and 
7.7 /zm relative to the CH out-of-plane bending mode at 
11.3 /im. Observed variations (shaded) and mean values 
(dashed) by Galliano et al. (2008). Full lines at various 
colors represent band ratios of PAHs with different hydro- 
genation, a(H/C), and excited by mono-energetic radia- 
tion with different hardness (E 7 , as labeled) and constant 
strength (flux, F). Band ratios derived from PAHs heated 
by the ISRF are indicated by filled circles and same color 
index of a. Beside other parameters the degree of ionisation 
is constant. 



Several teams who include PAH in MC codes (Wood 
et al. 2008; Baes et al. 2011, Robitaille 2011) use pre- 
computed emissivity profiles of the molecules, respectively 
P(T). The PAH libraries are computed by scaling a source 
spectrum of constant spectral shape (Draine & Li 2007; 
Compiegne et al. 2011) for which the interstellar radiation 
field of the Galaxy is often used (ISRF, Mathis et al., 1983). 
The libraries are therefore applicable only to radiation fields 
in which the strength is varied but its hardness remains 
close to the ISRF. In Fig. Q] we show that the band ratios 
of a PAH change within a factor of three when excited by 
photons of different hardness while they remain constant 
when excited by the ISRF. Therefore by including PAH 
into radiative transfer models the frequency of the photons 
exciting the molecules should be considered. This is in par- 
ticular true for hard photons of several eV well in excess of 
the mean photon energy of the ISRF (1.5 eV). In this case 
the pre-computed band strengths are incorrect by a factor 
close to what can be read from Fig. Q~] 

To cover most of the observed variations of the C-H 
and C=C band ratios as second parameter, the hydrogena- 
tion of the molecule given by a = [H] / [C] is used. We vary 
a between 0.05 and 0.4 which seems within the limit of 
realistic H coverage of such large molecules. We find that 
most of the observed values, with exception of those at 
-F(7.7)/-F(11.3) > 5, are reasonably covered without hav- 
ing changed the ionisation degree of the PAH (Fig. [1]). 

The reflection nebulae NGC1333, which is excited by an 
19,000 K photosphere, and NGC2068, which is heated by a 
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Fig. 2. Influence of stellar temperature and cut-off wave- 
length of the PAH absorption cross section on the IR emis- 
sion bands. Spectra are normalized to their far-IR peak. 

Table 1. PAH band ratios calculated for different stellar 
radiation fields and cut-off wavelengths of the PAH absorp- 
tion cross section, A cu t. 
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0.22 


0.66 
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75 


4.7 


0.34 


0.83 



(1) Blackbody temperature (hv < 13.6eV). 

(3) Cut-off wavelength in the PAH absorption cross section. 

(3) Average time between 2 photons absorptions (Eq.4, SK10). 

(4) Mean photon energy absorbed by a PAH (Eq.5, SK10). 

(5) Flux ratio at 6.2/^m and 11.3/^m. 

(6) Flux ratio at 7.7/zm and 11.3/im. 

11,000 K star, differ in the 7.7/11.3 /im band ratio by 30% 
(Uchida et ah, 2000). With the same stellar temperatures 
the ratios vary in the model by 15%, which can be increased 
varying the strength of the radiation field or the PAH pa- 
rameters. The cut-off in the PAH absorption cross-section 
is in the optical and moves only slightly to longer wave- 
length by increasing the size of the molecule. Laboratory 
measurements show that large neutral (Nq > 30), com- 
pact ions (Nc < 250) or nitrogenated PAHs may absorb up 
to the near-IR (Salama et al. 1996, Mattioda et al. 2005, 
2008). PAH emission spectra are shown in FigfS] using our 
ISM dust model. We vary the temperature of the star and 
A cu t- The stellar flux is constant at F = 100 ergs~ 1 cm -2 . 
Photons above 13.6eV are truncated. For A cut = 0.6^m 
(2eV) and in the 5,000 K stellar spectrum the PAH absorb 
every 140 s a photon with mean energy of 2.6 eV, whereas 
in the harder spectrum, the time between two absorption 



events is reduced to i a bs = 90 s and the mean photon en- 
ergy equals hv — 5.3 eV. The 7.7/11.3 band ratio changes 
by <^50% for both stellar temperatures (Table [T]), so that 
PAHs appear hotter in hotter stars. The picture changes 
when A cu t is increased to 1.8/xm (0.7eV). Then i a b s and 
hv changes as detailed in Table [1] and the band ratios vary 
with stellar temperature by 25% only. The choice of the 
PAH cross section has a strong influence on the emission 
spectrum. Unfortunately the cross sections vary by large 
factors from molecule to molecule. 



3. Radiative transfer in 3D including PAHs 

The radiative transfer problem is solved for arbitrary dust 
geometries by a MC scheme developed by Heymann & 
Siebenmorgen (2012). It includes grains large enough to be 
in thermal equilibrium with the environment. We first give 
a brief description of the procedure. Then we present an ac- 
curate method to include stochastic heated particles. The 
new code is verified against computations using ray-tracing 
methods. 



3.1. Large grains 

Our MC method considers advantages of different MC tech- 
niques (Lucy 1999; Bjorkman & Wood 2001). The compu- 
tational speed is increased by calculating the flight paths 
of many photons simultaneously. This is implemented by 
using vectorized hardware of either computer processing 
units (CPU) or graphical processing units (GPU). In the 
procedure a total of N — n ■ N v photon packages are emit- 
ted per second from a source of luminosity . In each of 
the N„ frequency bins there are n photon packages emitted. 
Each photon package has a constant energy e = L*/N. The 
model space is a three dimensional Cartesian grid (x, y, z) in 
which each cube can be divided into a number of subcubes. 
The flight path of each photon package is traced through 
the model. A package entering a cell interacts with the dust 
if the extinction optical depth along the path within the 
cell is At > —log(£), using unified random number £. The 
package is scattered if the dust albedo A > using ran- 
dom number otherwise it is absorbed. When the pack- 
age is scattered, it only changes direction. If the package 
is absorbed a new package of same energy but different 
frequency v 1 is emitted from the spot of absorption. The 
emission is isotropic. Each absorption event raises the en- 
ergy of the cell by e, and accordingly its dust temperature. 
In cells with extreme high optical depth photon packages 
may get trapped and for computational efficiency the mod- 
ified random walk procedure by Fleck & Canfield (1984) is 
implemented (Min et al. 2009; Robitaile 2010). The SED 
is computed by counting photon packets leaving the cloud 
towards the observer using a large beam. Images are com- 
puted by a ray-tracer where the radiative transfer of the 
line of sight from the observer, or a pixel of the detector 
plane, through the model cloud is solved using the MC com- 
puted dust absorption and scattering events. 
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3.2. Small grains 

Very small grains such as PAHs show strong fluctuations 
of their enthalpy state whenever they are hit by a photon. 
Their emission is generally treated statistically by intro- 
ducing a temperature distribution function P(T)dT which 
gives the chance of finding an arbitrary particle in a tem- 
perature interval [T, T + dT] . Absorption or emission leads 
to a transition from an initial Ui to a final U f enthalpy state 
of the system, which is represented by a transition matrix 
Au. For PAH we consider an analytic approximation of the 
enthalpy U (T) of graphite given by Chase et al. (1985). The 
enthalpy bins are of width AU{T) and binned in constant 
temperature intervals. Let for a given cell I, a number offti^ 
mono-chromatic photon packets of frequency v be absorbed 
by a PAH, then the transition matrix can be computed 



A fi = 



K V F 



ni tV hv 

: else 



if \U S - Ui - ni, v hv\ < \MJf 



(1) 



where K v denotes the PAH absorption cross section. In 
steady state the sum over all k-enthalpy states AjkPk = 
0, and P(T) is computed with the help of the rapid recur- 
sion formulae by Guhathakurta & Draine (1989). Examples 
of P(T) for mono-energetic excited PAHs are discussed 
by SK10. In some circumstances the distribution func- 
tion P(T) is rather fiat while in others it approaches a 5- 
function. Therefore we implement an iteration procedure in 
which the temperature grid of U(T) is carefully adapted to 
the specific case. This ensures that numerical inconsisten- 
cies in particular of the energy balance are avoided. Note 
that to reduce computing time such an iteration procedure 
is not included in the MC treatment by Bianchi (2008), who 
computes P(T) on a fixed temperature grid. 

Stochastic heated grains have been included in MC ra- 
diative transfer models with different levels of simplifica- 
tions. Kriigel (2008) describe two ways to treat stochastic 
heated particles in MC schemes. We follow the one where 
PAH are only excited by photons which are directly emitted 
or scattered from the source. A similar approach is taken 
by (Gordon et al. 2001; Misselt et al. 2001; Dullemond et 
al. 2007). A library approach is considered by Juvela & 
Padoan (2003), Wood et al. (2008), Baes et al. (2011), and 
Robitaille (2011), in which P(T) is pre-computed for some 
radiation fields such as a scaled version of the ISRF. This 
method reduces the high computational demand required 
to solve P(T) case by case, however at the expense of a 
possible computational error in particular for the derived 
PAH band strengths as discussed in Sect. [5J 

Previous MC codes which include PAHs do not include 
the physics of the dissociation process of the molecules; 
PAH photo-destruction by a single hard photon is neglected 
and soft X-rays are not considered. In our scheme we in- 
clude these effects as described below. Our code is optimised 
so that photons which are absorbed and re-emitted by dust 
are assumed not to heat the PAHs. This is a valid approx- 
imation as dust emits photons in the IR which are hardly 
absorbed by a PAH. This allows us to reduce storage and 
run time requirements. 



PAH are implemented in our code as follows: First, the 
source emits N photon packages and an equilibrium of the 
MC model is computed by considering the absorption but 
ignoring the emission of the PAHs. At this point only the 
number and frequency of the packets which are absorbed 
by a PAH are saved. Thereafter each cell with nf AH (zy) > 
becomes an independent source and we assume that it emits 
from its center the same number n of photons at frequency 
v' given by: 



e n {v)dv > £ / £ n {v)dv 



with random number £ and PAH emissivity 



PAH 



B„{T)P n (T)dT 



(2) 



(3) 



where K^ AH is the PAH absorption cross section. 

One extra complication occurs for cells where PAH dis- 
sociate. For nanometer-sized grains such as PAH we con- 
sider photo-destruction. It is estimated that a PAH with 
Nc carbon atoms is destroyed (Omont, 1986), when the 
absorbed energy AE exceeds a critical energy E CI i t : 



AE > £ cri t = ^ (eV) 



(4) 



The energy AE, may be either due to absorption of i) many 
soft photons or ii) whenever the molecule is hit by a single 
hard photon. The first condition is only relevant at short 
distances from the source while the second acts irrespec- 
tive of that distance. If EqfJ] is fulfilled and PAHs of cell i 
dissociate the nf AH (z/) absorbed photon packets are again 
emitted from the star. We further simplify the method and 
ignore PAH absorption during the flight path of these pho- 
tons through the cloud. 



3.3. Verification 

For equilibrium heated grains the MC code is veri- 
fied against various one and two dimensional bench- 
marks (Heymann & Siebenmorgen, 2012). For stochas- 
tic heated particles we compare results computed by MC 
against the spheric symmetrical radiative transfer codes 
by Siebenmorgen (1993) and Kriigel (2008). In the latter 
the radiative transfer problem is solved by a ray tracing 
method. We find that the SED computed by the MC pro- 
cedure and the ray tracing code agrees typically to within 
a few percent (<5%). One such comparison is shown in 
Fig. [3] The dust sphere is heated by a star with luminosity 
of 1L Q and photospheric temperature of 6,000 K, between 
r in = 3AU and r out = 90 AU a constant dust density of 
1.85 x 10 -19 g/cm 3 is used. For the large grains we take 
ISM type dust with parameters as in SHI 2 and PAHs as of 
Sectj2] This gives a dust extinction cross section in the V 
band of K^ xt = 43, 000 [cm 2 /g-dust]. The dust mass of the 
cloud is 9.53 x 10 _7 M© and the optical depth between r in 
and r out is ry = 10. The overall SED, including the far-IR 
emission and the PAH feature strengths of both radiative 
transfer computations are in good agreement (Fig. [3]). 
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Parameter 


T Tauri 


Herbig Ae 


Stellar luminosity L, [Lq] 
Stellar mass M, [M ] 
Photosphcric t6nip6r3,tur6 T* [K] 

Column density E(r) - T ^ AU) ( a ^)t [cm 2 /g-dust] 

Vertical optical depth rx(lAU) 

Dust density in halo Phaio [g-dust/cm 3 ] 

Distance D [pc] 

Inner disk radius ri n 

Outer disk radius r out [AU] 


2 
1 

4,000 
r < 1 AU 
r > 1 AU 
1C 

or 1. 

evap 

22.5 


50 
2.5 

10 000 

: 7 = 0.5 
: 7 = -l 
,000 

5 x 10" 18 
50 

oration 
40 



10.0 F 




wavelength (|im) 

Fig. 3. Comparison of the SED of a stellar heated, optically 
thick (tv = 10) dust sphere of constant density computed 
with a ray tracing method (Kriigel 2008, full line) and the 
MC code of this work (dashed) . The flux ratio between both 
methods is displayed in gray. 

4. Protoplanetary disk models 

Protoplanetary disk models are set up similarly to those 
described in SHI 2. We consider both photospheric heating 
and hard radiation components of the star. In addition to 
large grains we include PAHs. First, the set up of the hy- 
drostatic disk models are described. Then our treatment 
of the PAH destruction in the disks is outlined and results 
compared for disks computed in ID for flaring disks and in 
2D by the MC scheme. Models are applied to disks heated 
by a T Tauri star and a Herbig Ae star. Halos above the 
disks are considered. The surface brightness in the PAH 
bands is discussed. 

4.1. MC disk models 

As fiducial model of the T Tau star we take a mass of 1M Q , 
a luminosity of L* = 2L Q , and a photospheric temperature 
of 4,000 K. For Herbig Ae star we use 2.5M , 5OL and 
10 4 K (van den Ancker et al. 1997). Soft X-rays (50 eV - 



2keV) are included with luminosities up to 10~ 4 L*. The 
distance is set to 50 pc. 

For large grains we take the fluffy grain model of SH12. 
It includes a power-law size distribution: n(a) oc a -3 5 , 
with particle radii between 320A < a < 33/im of com- 
posite grains with volume fraction of 34% silicates, 16% 
carbon, and 50% vacuum. Dust abundances (ppm) are 31 
[Si]/[H] and 200 [C]/[H] (Asplund et al., 2009). Optical 
constants are by Draine (2003) for silicates and Zubko 
et al. (1996) for carbon. Absorption and scattering cross- 
sections and the scattering asymmetry factor is computed 
with the Bruggeman rule. This gives a total mass extinction 
cross section for the large grains in the optical (0.55/j,m) of 
j^ext _ 4 ; 000 [cm 2 /g-dust] and a gas-to-dust mass ratio 
of 130. PAH parameters are taken from Sect [51 Kinetic en- 
ergy losses by e _ are important when dust interacts with 
energetic photons of critical energy >100eV. This effect is 
not included in the Mie theory and will result in an extra 
~ 1/v decline of the absorption cross section; for detailed 
computations and size dependence of the critical energy 
see Dwek & Smith (1996). This correction is applied to the 
cross sections of all particles. 

The disk extends inwards up to the point where porous 
grains reach 1,000 K, or equivalently small interstellar 
grains would be at about 1,500 K. PAH destruction is com- 
puted as discussed in Sect. 14.21 The outer disk radius is 
r out = 22.5 AU for T Tau and 40 AU for Herbig Ae disks. 
Further out the midplane temperature drop below 30 K 
and other heating mechanisms may become important. The 
PAH emission is well constrained within a fraction of r out . 
Initially we assume that the disk is isothermal in the verti- 
cal direction z and the density is given by 

P(r,*) = y|§e-«V»*" (5) 

with scale height H 2 = kT m idr 3 /GM*m, surface density 
S(r), molecular mass m — 2.3m p and midplane temper- 
ature, T m jd, for which we use a power-law as an initial 
guess. The height of the disk is set initially to zq = 4.5H. 
In models with pure disks the density phaio above z > z Q is 
zero and constant when a halo is considered (Tabled]). The 
surface density is adjusted to the optical depth in vertical 
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direction: t± = Y,{r)Ky = tiau (r/AU) 7 with 7 = — 1 for 
r > IAU and 7 = 0.5 otherwise (Min et al., 2011). At 1 AU 
the vertical optical depth from the surface to the midplane 
is tj_(1AU) = 10,000. This translates to a surface density 
of S(1AU) = 5g-dust/cm 2 for disks without PAHs, which 
is close to estimates for the minimum mass of the early so- 
lar nebulae (Hayashi, 1981), and S drops by a factor 4.4 
when PAHs are considered. 

The disks are in hydrostatic equilibrium in the verti- 
cal direction so that gravitational force is balanced by the 
pressure gradient 

z GM* _ 1 dP 
r r 2 p dz 

with pressure P — pkT(z)jm. We solve Eq. [5] and Eq. [B]by 
an iterative scheme: After a first set up of the disk struc- 
ture p(r, z) by Eq.[5l the temperatures T(r, z) are computed 
with the MC code. By inserting the temperatures into Eq.[5] 
a new density structure p(r, z) is found. Then with a new 
MC run T{r, z) are updated and the procedure repeated. 
After less than a dozen of iterations a stable disk configu- 
ration is found. 

4.2. PAHs destruction and X-rays 

Young low-mass stars, and to a lesser degree intermediate- 
mass stars, are X-ray emitters. Here we study the influ- 
ence of hard photons on the appearance of PAH bands 
in disks. For the destruction of PAHs the process known 
as unimolecular dissociation is applied (Tielens, 2005). A 
simplified scheme which allows estimating of the distance 
from the star where the molecules become photo-stable is 
described in SK10. The threshold energy AE (Eq. g|) to 
dissociate PAHs by multi-photon absorptions is given by: 

AE = N c J^^du >E clit (7) 

where L v is the frequency dependent luminosity of the star, 
k v is the absorption cross section per C atom of a PAH 
and r„ is the optical depth from the star to the spot of ab- 
sorption in the cell. Hard photons can dissociate PAHs at 
all distances from the source in a time short compared to 
the life time of the disk. Therefore a PAH survival channel 
is introduced for which, as a result of turbulence, vertical 
mixing of the dust within the disk is assumed. From Spitzer 
spectroscopy Keller et al. (20008) suggest that PAH may be 
formed by evaporation of icy grain mantles and constantly 
destroyed and replenished. We crudely simplify the picture 
and assume that for the size of the largest eddies the aver- 
age turbulent velocity grows linearly with the sound speed. 
The turbulent velocity is parametrized vj_ = \fa c s with 
sound speed c s . Estimated values for a are often used but 
highly uncertain. They range from a ~ 0.0001 up to 0.1 
(Schrapler & Henning 2004, Youdin k Lithwick 2007), if 
not otherwise mentioned we use a ~ 10~ 3 . The influence of 
the choice of a on the PAH survival is discussed in SK10. 
For a more realistic physical description of the turbulent 
motions the MC treatment has to be combined with three 



dimensional hydrodynamical calculations (Armitage 2011, 
for a review). Excitation of PAHs by stellar photons is to 
good order within the extinction layer of optical radiation 
and has thickness I ~ H/2. For vertical velocity v± — £/t exp 
a PAH is exposed to radiation by a time i CX p. Using the 
destruction time ides, as the time between two hard pho- 
ton absorption events, the PAH will survive if ides > iexp- 
Translated into a velocity condition the PAH hit by photons 
with -E cr it = hv CT n will be destroyed if (Eq.23 in SK10): 

00 

v± < v CT = £/t dcs = — — / — dv (8) 

fcrit 

with ft = l/H ~ 0.5. Both PAH destruction conditions 
of Eq. [7] and Eq. [8] are straight forward to implement into 
radiative transfer codes which are based on ray-tracing. 

By implementing the hard radiation components into 
the MC code we enter into numerical challenges. The 
equidistant energy binning of the photospheric component 
is achieved with 2,000 frequencies. However, considering 
stars with non-negligible far-UV, extreme-UV and X-ray 
components means increasing the frequency grid to 250,000 
bins at the expense of a run time increase by a similar fac- 
tor (>100). Therefore we need to simplify the treatment 
of PAH dissociation by hard photons in the MC code. The 
energy balance of the dust disk is dominated by the pho- 
tosphere and the hard radiation components are of small 
enough luminosity that they play a minor role in the over- 
all disk structure and the temperature of the midplane. 
Hard radiation plays a major role in computing the PAH 
destruction zone. The far-UV photons have a mean en- 
ergy below -E cr it and do not contribute to the condition of 
Eq. [7] Most of the extreme-UV photons will be absorbed 
by the gas located well above the height of dust extinction 
layer (Ercolano et al., 2010). Therefore it is safe to ignore 
far-UV and extreme-UV radiation components. The de- 
tailed spectral variation of X-rays is not relevant for the 
destruction of PAHs. This allows us to represent X-rays by 
a mono-chromatic energy which we choose to be at 1 keV. 
This choice ensures that each X-ray absorption event leads 
to complete destruction of the molecule. To apply the PAH 
destruction condition of Eq|S]for each cell only the total op- 
tical depth, rikeV remains to be computed. Therefore the 
survival of PAHs in a cell can be estimated before the time 
consuming MC run. In cells where PAHs are destroyed, by 
the condition of Eq. [5J there will be no photon absorption 
by PAHs. 

4.3. T Tauri disks 

We compute the SED of hydrostatic disks with PAHs and 
hard radiation components for T Tauri stars applying three 
different disk geometries. First, we consider a flaring disk 
in which the radiative transfer is solved only in the vertical 
direction. Then we apply the MC code to disks assuming- 
axial symmetry. In the MC models the region above the disk 
is either free of dust or filled in a halo like configuration. 
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Fig. 4. Mid IR emission of T Tau disks with X-ray lumi- 
nosity of Lx = 10~ A L* (dashed), Lx = 10~ 6 L„» (dashed), 
and Lx = (full line) . SEDs are computed for flaring disks 
in 'ID' (top) and with MC (bottom). 



We apply the one dimensional (ID) radiative transfer 
code developed by Kriigel (2008) of a hydrostatic and ge- 
ometrically thin disk. The program is applied to the emis- 
sion of PAHs from disks of Herbig stars by Habart et al. 
(2004) . The disk is symmetric with respect to the midplane 
at z = 0. The density structure is given in cylindrical coor- 
dinates p(r, z). Light from the star falls on the disk under a 
grazing angle a gr so that from everywhere on the disk sur- 
face the star is visible (Armitage, 2007). For the ID models 
we use the flaring geometry by Chiang & Goldreich (1997) 
set to 3° < a gr oc r 2 / 7 < 7°. We present results where the 
ID disks are heated by four radiation components of the 
star (SK10): the photosphere with luminosity of 0.99 L* and 
black body temperature of 4,000 K, far-UV with 1%L* at 
15,000K, extreme-UV with 0.1% L* at 300,000K and X- 
rays of up to 2keV with Lx = 10~ 36 L* and a spectrum 
oc v 2 . We implement into the code the PAH destruction 
mechanism following conditions of Eq. [7] and Eq. [SJ 

The SEDs of T Tauri disks including PAHs and heated 
by hard radiation are computed also by the MC scheme 



Fig. 5. IR emission of T Tau disks with and without dust 
envelope, viewed nearly (~ 20°) edge-on. Heating includes 
X-ray luminosities as indicated. 

for which parameters of Table [5] are taken. We find that 
the appearance and non-appearance of PAH bands is pre- 
dicted in a similar way in both geometries provided by the 
ID and the MC model. As shown in Fig. [H if one consid- 
ers only the photosphere strong PAH bands appear in the 
spectra. They weaken when including hard radiation fields 
of Lx = 10 -6 L* and PAH do not contribute to the emission 
spectrum for Lx > 10 -4 L*. The latter X-ray luminosity 
can be considered as a lower limit for most T Tauri stars 
(Preibisch et al., 2006). Our computations corroborate ear- 
lier estimates that the non-appearance of PAH bands in 
spectra of T Tauri stars is connected to the destruction of 
the molecules by hard radiation (SK10). We verified that 
X-ray luminosities have little impact on the midplane tem- 
perature, nor do they impact the detailed surface struc- 
ture of the disks. The ring-like structures, their shadows 
and gaps in the disk surface are preserved as computed by 
SH12. 

4.4. Halos 

Residual dust may be left or replenished in a halo at high 
altitudes above the disk (Vinkovic et al. 2006; Krijt & 
Dominik 2011). Early models added spherical distributed 
and optical thin dust emission on top of SEDs computed 
for pure disks (Miroshnichenko et al., 1999). The MC code 
allows a self-consistent treatment of the radiative transfer 
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Fig. 7. Mid-IR images at 11.3/im of a disk with (middle) and without (right) PAHs heated by a Herbig Ae star and 
viewed at 30° from the midplane. Images are normalized to the stellar flux. Color bars are in log(Jy/arcsec 2 ). Left: 
Surface brightness distribution of the images measured as cut from the star. 



in such a geometry. The halo is filled by large grains and 
PAHs and is implemented into the code assuming a dust 
density of p(r,z) > phaio = 1-5 x 10~ 18 (g-dust /cm 3 ) . This 
gives an optical depth of the halo of Ty ~ 1.5. Resulting 
SEDs are shown in Fig. [5l The X-ray luminosity is varied 
using Lx = 0, 10 _6 L*, and 10 _4 -L*. The dependence of 
the PAH band strengths on the X-ray luminosity is similar 
to models without halos (Fig. [4]) . The halo models appear 
warmer and have a factor ~ 3 stronger IR peak emission 
and at <6/im less near IR flux than pure disks. The halo 
provides additional heating of the disk on top of the direct 
stellar light and therefore their disk midplane temperature 
is warmer and shows less structure. 

4.5. Herbig Ae disks 

We present disk models for Herbig Ae stars with parame- 
ters of Table [2] Resulting SEDs are computed by the MC 
scheme and are shown for face-on and edge-on views in 
Fig. El The optical depth of the halo is Ty ~ 1.5. The influ- 
ence of hard radiation on the PAH destruction is studied 
by applying a X-ray luminosity of Lx = and 10 L*, 
which exceeds the range of 10 -5 ... 10~ 7 L* typically found 
in Herbig Ae stars (Stelzer et al., 2006). Irrespective of these 
X-ray luminosities PAH bands are visible in the Herbig Ae 
disks. This is contrary of what we found for T Tauri disks 
(Sect. l4~3j) . At the inner rim of the disks the absolute X-ray 
flux is more than a factor 10 lower in Herbig stars than in T 
Tauri stars and PAH survive. If one lowers a by one order 
of magnitude the PAH band-to-continuum ratio is reduced 
by 30% and for a = the PAH signatures are absent. We 
verified that the up and downs in the midplane temperature 
as well as the ring-like structures, the shadows and gaps in 
the disk surface are preserved as computed by SH12, irre- 
spectively of the X-ray luminosity or the strengths of the 
PAH emission. 

When the disks are viewed face-on, the PAH feature- 
to-continuum ratio becomes smaller and even more so when 



a halo is present. In the latter configuration the emission 
by large grains outshine the PAH contribution. We noticed 
a similar effect in disks of T Tauri stars in which the X-ray 
luminosities are small enough for PAH survival. Therefore, 
the geometrical distribution of the dust effects the detection 
probability of PAH features in disks. 

The stochastic heating of PAH is rather independent of 
the distance to the star, a fact first noticed in reflection 
nebulae (Sellgren et ah, 1985). We study how this effect in- 
fluences the appearance of pure Herbig Ae disks in the mid- 
IR. The surface brightness distribution in the 11.3/im band 
is displayed in Fig.[7]for disks with and without PAHs. The 
images with PAH emission are brighter than those without 
PAHs at distances between 6AU and 10 AU. Otherwise, 
the emission structure remains the same, except for some 
broadening of the inner gap in the disk, outside the dust 
evaporation zone. In both cases the dominant contribution 
of the emission is coming from the central 2-3 AU. Due 
to the brighter optical luminosities of Herbig Ae stars the 
PAH are excited at large distances from the star. There 
the turbulent velocity is larger than the critical velocity for 
PAH destruction by X-rays, so that the molecules survive 
(Eq. [5]) . For a close to the surface brightness distribution 
approaches that of disks without PAHs. 

5. Conclusion 

We study the dependencies of the PAH band ratios by 
changing the hydrogenation coverage of the molecules and 
the spectral shape of the exciting radiation field. PAH band 
ratios vary by a factor three when excited by photons of 
different hardness. They remain almost constant when ex- 
cited by the ISRF. Therefore by including PAH into radia- 
tive transfer models the frequency of the absorbed photons 
should be considered. 

MC radiative transfer models which treat stochastic 
heated grains implement different levels of simplifications: 
Temperature distribution functions are often used which 
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Fig. 6. IR emission of a Herbig Ae disks with (magenta) 
and without (black) dust envelope, viewed nearly face- 
on (top) or edge-on (bottom). Results for different X- 
luminosities as indicated. 



are pre-computed for a constant hardness of the exciting 
radiation, photo-dissociation of the PAHs is not included 
and hard radiation components are ignored. The MC code 
presented in this paper considers heating of the PAH by in- 
dividual photons and photo-destruction of the molecules. 
We compare our PAH treatment in the MC scheme against 
a radiative transfer program which is based on ray-tracing. 
The PAH emission computed by both codes is in good 
agreement. 

The PAH emission in protoplanetary disks is studied for 
T Tauri and Herbig Ae stars. Dust in the disk is configured 
either in a slab geometry, in which the radiation transport is 
solved in vertical direction (ID), or by the MC procedure 
assuming hydrostatic and radiation balance (SH12). The 
disks are heated by the stellar photosphere and soft X-rays. 
Hard photons dissociate PAHs in the disk within a time 
that is short compared to the lifetime of the disk. Therefore, 
following SK10, a PAH survival channel is introduced in 
the radiative transfer programs for which, as a result of 
turbulence, vertical mixing of the dust within the disk is 



assumed. We find that the appearance and non-appearance 
of PAH bands is predicted in a similar way by the ID and 
the MC disk models. 

We find that stars with X-ray luminosities > 10 4 L do 
not display PAH emission features in disks around T Tauri 
stars whereas PAH are present in disk models of Herbig Ae 
stars because their absolute X-ray flux is more than a factor 
10 lower. Due to the brighter optical luminosities of Herbig 
Ae stars, the PAH are excited at large distances from the 
star. There the turbulent velocity is larger than the critical 
velocity for PAH destruction, so that the molecules sur- 
vive and PAH emission can be detected. The dependence 
of the X-ray luminosity on the PAH band ratios is similar 
in disk with and without halos. However, in disks with ha- 
los the PAH band ratios are smaller than in disks without 
halos. The mid-IR surface brightness appears brighter in 
the planet forming region of disks with PAH emission than 
in those without PAHs. 

Acknowledgements. We are grateful to Endrik Kriigel for helpful com- 
ments and for providing his ID radiative transfer codes. 
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